Multiscale estimation of GPS velocity fields

S U M M A R Y We present a spherical wavelet-based multiscale approach for estimating a spatial velocity field on the sphere from a set of irregularly spaced geodetic displacement observations. Because the adopted spherical wavelets are analytically differentiable, spatial gradient tensor quantities such as dilatation rate, strain rate and rotation rate can be directly computed using the same coefficients. In a series of synthetic and real examples, we illustrate the benefit of the multiscale approach, in particular, the inherent ability of the method to localize a given deformation field in space and scale as well as to detect outliers in the set of observations. This approach has the added benefit of being able to locally match the smallest resolved process to the local spatial density of observations, thereby both maximizing the amount of derived information while also allowing the comparison of derived quantities at the same scale but in different regions. We also consider the vertical component of the velocity field in our synthetic and real examples, showing that in some cases the spatial gradients of the vertical velocity field may constitute a significant part of the deformation. This formulation may be easily applied either regionally or globally and is ideally suited as the spatial parametrization used in any automatic time-dependent geodetic transient detector.